Ga-NMR local susceptibility of the kagome-based magnet 
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Abstract 

We report a high-T Ga-NMR study in the kagome-based antiferromagnetic 



^ . compound SrCr9pGai2-9pOi9 (.81 < p < .96), and present a refined mean- 

o, 

, ^,1 field analysis of the local NMR susceptibility of Cr frustrated moments. We 



find that the intralayer kagome coupling is J = 86(6) K, and the interlayer 
coupling trough non-kagome Cr moments is J' = 69(7) K. The J' /J = 0.80(1) 
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The kagome-based compound SrCr9pGai2-9pOi9[SCGO(p),p < 1] has been receiving 
considerable attention as a model system for geometrically frustrated physics [0. Mag- 
netic properties of SCGO(p) arise from Heisenberg Cr'^"'" ions (5 = |). Despite macroscopic 
susceptibility (xmacro) measurements which indicate that a strong antiferromagnetic inter- 
action couples neighboring spins {Qmacro ~ 500 — 600 K), spin-spin correlation never exceeds 
twice the Cr— Cr distance 0. No classical Neel long range order occurs. On the contrary, 
SCGO(p) continues to display Curie- Weiss behavior well below Qmacro as suggested by the 
linear T-dependence of X^acroy extending mean-field predictions to an unusual temperature 
domain (T << Qmacro)- 

Geometric frustration arises from interacting spins on kagome networks. However, the 
crystal structure is not pure kagome since Cr^+ ions occupy both kagome [Cr(12/i;)] and non- 
kagome [Cr(2a),Cr(4/„j)] sites, as shown in the inset of Fig.|l]. Cr(2a) sites link two kagome 
layers. The kagome-Cr(2a)-kagome slabs are isolated from one another by Cr(4/t,j)— Cr(4/^j) 
pairs which Lee et al. showed to form spin singlets with a binding energy of A = 216 K 
0]. They also suggested that the kagome-Cr(2a)-kagome slab may be regarded as a quasi- 
2D highly frustrated system, which they named pyrochlore slab. However there is a lack of 
experimental evidence concerning couplings which ensure this slab structure. Coupling of 
the two kagome planes through Cr(2a) is not known, nor is the intralayer kagome coupling. 
One way to shed light into this problem, would be a refined high-T analysis of the magnetic 
susceptibility of Cr frustrated moments. Unfortunately this is prevented in Xmacro by the 
fact that there are contributions from all magnetic sites, frustrated and non-frustrated. 
Moreover, Xmacro is also sensitive to disorder as non-magnetic Ga^^ ions are always present 
on Cr^+ sites (the p = 1 compound is not stable) Q. As we show here, such an analysis 
can be done by ^^Ga and ^^Ga (/ = |) NMR since Ga nuclei on site 4/ directly probe both 
Cr(12A;) and Cr(2a) moments (see. inset Fig.|l]). By varying Ga/Cr substitution from 4% to 
19% {p = .81, .89, .96), we are able to separate and identify on a firm experimental ground 
the substitution-related susceptibility from the intrinsic frustrated susceptibility of the slab 
I^J^. We present here high-T results and mean-field analysis of the slab susceptibility, 



which enables us to evaluate J' and J, respectively the Cr(2a)— Cr(12A;) coupling and the 
Cr(12A;)-Cr(12fc) coupling. We find J' = 69(7) K and J = 86(6) K. The J'/ J = 0.80(1) 
ratio confirms the common belief that the frustrated entity is a pyrochlore slab. 

All samples are ceramics and were synthesized by solid state reaction as detailed elsewhere 
1^. Reaction products were checked by standard x-ray diffraction and refined by neutron 
diffraction (see Table I). Xmacro yielded results in agreement with previous concentrations 
reported in literature {Qmacro = 439(22), 501(14), 608(14) K for p = .81, .89, .96 respectively) 

i- 

An extended NMR analysis of the spectrum was presented in previous papers [^^f. 
Three different Ga NMR sites can be resolved: Ga(4e), Ga(4/) and Ga substituted on Cr 
sites [Ga.{sub.)]. Ga(4/) is at the heart of frustrated physics (inset Fig.|l|). It is exclusively 
coupled to Cr(12A;) and Cr(2a) by a Ga— O— Cr hyperfine interaction, and is the object of 
the present study. 

A set of spectra ranging from 150 K to 410 K were performed in a ~7 Tesla permanent 
magnet. We used a 7r/2, vr sequence with spacing r = 18 /is between pulses, and as the 
applied frequency was swept over a ~1 MHz range, with a step interval of 75 kHz, spin 
echo signals were recorded. The overall spectrum was then obtained by summing all the 
spin echo Fourier transforms [^. Another set of spectra ranging from 15 K to 200 K were 
obtained in a sweep field setup, i.e. by varying an external magnetic field at a constant 
applied frequency 40.454 MHz, and recording the integral of the spin echo signal. Again we 
used a conventional 7r/2, vr sequence, but with spacing r = 28 fis. 

A typical set of high-T (T > 100 K) field sweep spectra for the p = .96 sample, recorded 
around 3 Tesla, is reported in Fig.||. In this field window two lines are unambiguously 
resolved. Ga(4/) is the main line and the hump at its right is Ga.{sub.). The negligible 
contribution of Ga(sub.) to the spectrum refiects the nearly stoichiometric value of the 
sample. The sharp non-symmetric features of the Ga(4/) line originate from T-independent 
quadrupole effects (^^uqlAf) ~ 2.9 MHz, 77 ~ 0) related to local charge environment around 
the nucleus. Upon cooling the NMR line shifts to the left with respect to its zero reference 



position, without any appreciable broadening. The T-dependence of the shift (K) stems 
from the Ga(4/) magnetic environment. K provides a direct measurement of the magnetic 
susceptibihty of the Cr(12/i:) and Cr(2a) ions (to be detailed below). 

The temperature dependence of K~^ is presented in Fig.^ for all samples. The high-T 
data (T > 100 K) follows a Curie- Weiss law as suggested by the linear variation of K^^. By 
extrapolating K'^ to we extract the NMR Curie- Weiss temperature {Qna4r)- We evaluate 
it to be 453(50) K, 469(44) K, and 484(40) K for p = .84, .89, .96 respectively, of the same 
order of Qmacro, but somewhat lower. The increase in Qnmr with p reflects a better lattice 
coverage as one would expect from mean-field theory. Here again Curie- Weiss behavior 
extends well below Qnmr- Further lowering T, K in all samples eventually deviates from 
Curie- Weiss behavior, fiattening and even decreasing below T = 40 — 50 K. We address the 
reader to other publications for a detailed discussion on the low-T behavior of K p,p . 

We now turn to discuss the high-T behavior of K within a mean-field analysis. First we 
would like to point out that Cr(12A;) and Cr(2a) do not have same magnetic environment. 
Cr(2a) has 6 nearest neighbors of Cr(12A;), whereas Cr(12A;) has 5 nearest neighbors (4 
Cr(12A;) and 1 Cr(2a)). Furthermore Cr(12A;)-Cr(12A;) coupling (J) and Cr(12A;)-Cr(2a) 
coupling (J') are a priori not the same. In these regards our mean-field analysis is similar 
to that performed in ferrimagnetic compounds which bear inequivalent magnetic sites. 

To justify this analysis we must argue how J and J' are the only relevant couplings in the 
slab. A study of a series of chromium oxides indicates that the Cr-Cr interaction is dominated 



by direct exchange |jTO|. It was shown that such an interaction varies rapidly with Cr— Cr 
distance (and Cr— O— Cr angle), being ^ when dcr-Cr ^ 3.1 A. Therefore only interactions 
between Cr(12A;)-Cr(12A;) {dcr-Cr = 2.895 A) and Cr(12A;)-Cr(2a) {dcr-Cr = 2.971 A) are 
relevant in the frustrated unit, with < J' < J. Cr(2a)— Cr(2a) interaction is ~ as 
dcr-Cr = 5.80 A. Note that the kagome layer is distorted, therefore J represents the average 
coupling between Cr(12A;). 

On the basis of these considerations, the frustrated spin Hamiltonian is 



H = J 2_^ (Si2k)i ■ (Si2k)j + £(Si2k)i ■ C^2a)j 

(M) 

where we have introduced < e = J'/ J < 1. The first term is the isotropic Heisenberg 
interaction between spins 012k on the kagome layer. The second term is the interaction 
between ouk and o 2a spins. The summation is taken over neighboring spins of the slab. 

Using a mean- field approach the Xi2k and X2a susceptibilities per site of Cr(12A;) and 
Cr(2a) can be written 

^ /^e// . P2aS{S + l)eJ 

^''"^ 3kBTf{Ty 3kBT '' ^^ 

and 

^ /^e// . 6pukSiS+l){2/3-E)J 

where /(T) = l + 4:pi2kS{S + l)J/3kBT — 6pi2kP2a{£JS{S + l)/3kBTy- The two expressions 
are different and reflect the fact that the two sites have different magnetic environments. 
The effective magnetic moment is fieff = 3.83fiB, typical of S* = 3/2 Cr'^+ ions |ni|. pi2k,P2a 
are the occupation of the Cr(12A;) and Cr(2a) sites known from neutron refinement (cf. Table 
I). The only unknown parameters are J and e. In Fig.^ we show the evolution of Xi2k and 
X2a when increasing e, with a fixed value J. In the limit e = 0, Xi2k follows a Curie- Weiss law 
with a kagome Curie- Weiss temperature 4:pi2kJS{S + l)/3. Whereas X2a follows a Curie law 
as we have freed the Cr(2a) spin from its magnetic environment. As e is increased, the pure 
paramagnetic behavior of X2a progressively dies out, while Xi2k is little affected. Of interest 
is when the kagome-Cr(2a)-kagome slab behaves as a pyrochlore slab. By pyrochlore slab 
we mean that both magnetic sites are identical Xi2k = X2a, i-e. e = Eq = 4:pi2k/{Gpi2k —P2a)- 
For an homogeneous substitution {pi2k = P2a = p) ^o = 0.8, reflecting the fact that the 
effective exchange field on Cr(2a) (6J') and on Cr(12A;) (4J + J') are then identical. We 
stress that this relation is independent on distortion on the kagome layer, as it only involves 
the average exchange constant J. 

We may express the NMR shift K using the mean-field expressions (|ID and (|]) 



K OC {9pi2kXl2k + 3Sp2aX2a) (3) 

where we have introduced 6 = A' /A. A, A' are the hyperfine constants which are respectively 
a measurement of the Cr(12/c),Cr(2a) moment at the nuclear site, and just influence the 
thermal behavior of K by their ratio 6. The two terms of expression (^ are multiplied by 
the number of neighboring Cr(12A;) and Cr(2a) linked by the oxygen ions to the Ga(4/) 
nucleus, respectively 9 and 3. We can fix a lower bond to 6 by some simple structural 
remarks. The hyperfine bridge Ga(4/)— O— Cr(12/c) yields bond lengths oi dca{4.f)-o = 1-866 
A and do-cr{i2k) = 2.052 A. The hyperfine bridge Ga(4/)— O— Cr(2a) yields bond lengths 
of dGa{4:f)-o = 1-866 A and do-cr{2a) = 1-971 A. The small difference in length and the 
similarity for angles between these bonds leads to a stronger coupling to Cr(2a) then to 
Cr(12A;), and 5 > 1. A more precise calculation suggests 5 ~ 2 0. It is instructive to 
compare K to the expression of the susceptibility of the kagome-Cr(2a)-kagome slab that 
one would obtain from a macroscopic measurement 

Xslab = 6pi2fcXl2fc + P2aX2a 

where Xi2k and X2a are now weighted by 6 and 1, respectively the number of Cr(12fc) and 
Cr(2a) per formula unit of SCGO. K is an independent measurement of the slab suscepti- 
bility, but with a different weighting than Xsiab- Using 5 ~ 2, we see that X2a, Xi2k impact 
i^ in a ratio of 6/9 instead of 1/6 for Xsiab- K on the contrary of Xsiab is very sensitive to 
temperature dependence of X2a, meaning that it is also very sensitive to e. 

The shift K reseals three unknown parameters J, e, and 5. Our aim is to determine 
their values by fitting for T > 100 K the shift-data presented in Fig.^]. There are an infinity 
of (e, J, 5)-set of values which reproduce the linear variation of K^^ with T. In Fig.^ we 
present the set of (e, J, 5)-values obtained when fitting K~^ of the p = .96 sample with 
expression (^. As shown, for a given value of e (abscissa) there is an associated value 6 
(upper panel) and J (bottom panel). In order to seize the evolution of this set of parameters 
it is instructive to examine the dependence of 6 on e. In the limit e — * 0, we see that 
5 ^ 0. In fact when e = 0, X2a follows a Curie law. Therefore to reproduce the linear 

6 



variation of K~^ there is only one possibility: uncouple the Ga(4/) nuclei from Cr(2a), 
setting (5 = 0. As e is increased, 6 monotonically increases since Curie contribution from 
X2a dies out. But when e = Eq expression (^ reduces to i^ oc {9pi2k + 35p2a)xi2fc, since in 
this case Xi2k = Xia = Xsiab/{Qpi2k + P2a) ■ The T-dependence of the shift is then completely 
determined by the {e, J) pair of values, whatever the value of 6. When e > eo,we recover the 
monotonic increase in S. Of interest is also the slow monotonic decrease of J with increasing 
e. Since < J' < J, i.e. < e < 1, we conclude that 80 K< J < 100 K. However only the 
set of values for which 6 > 1 bear a physical solution. Two possibilities subsist: (i) the set 
of values close to e ^ Eq (ii) the set of values for e > A. Case (ii) is a non-physical solution. 
This leads us to one of our major findings: the high-T linear behavior of K~^ reflects the 
physics of a pyrochlore slab, with the consequence that K oc Xsiab- 

Results for all samples are presented in Table I. Errors are governed by error bars on pi2k 
and p2a from neutron diffraction data. Coupling constants are found independent on lattice 
coverage given the negligeable variation of lattice parameters with substitution. Turning 
to the value of J, it is very satisfying to find a coupling constant close to the one ob- 
served in Cr203, where Cr'^'^ ions have a local octahedral environment {dcr-Cr = 2.890 A, 
Cr-0-Cr= 93.1°) similar to Cr(12A;) {dcr-cr = 2.895 A, Cr-0-Cr= 93.8°). There, it 



was established that the exchange coupling between neighboring Cr is 77(3) K ||T2[. We 
also performed a consistency test by calculating the expected values of Qnmr using re- 
sults reported in Table I. The analytical expression of Qnmr is obtained from the high-T 
limit (T » S{S + 1)J) of expression (§. We find 456(56), 460(33) and 495(35) K for 
p = .81, .89 and .96 respectively, in agreement with experimental values of Qnmr- Finally 
by using the susceptibility of Cr(4/i,j)— Cr(4/i,j) pairs ||^, we are able to reproduce the high-T 
Xmacro— data where the susceptibility Xsiab adds on to the Cr(4/^,j) and to the defect induced 
susceptibilities. 

In conclusion, we have evidenced that the Curie-Weiss behavior of the frustrated slab 
susceptibility extends to much lower temperatures than the average Cr— Cr interaction. 
A high-T mean-field analysis allowed us to determine couplings involved in the frustrated 



physics. We confirm the common behef that frustration in SCGO(p) arises from a pyrochlore 
slab. 

We wish to acknowledge fruitful discussions with A. Keren, C Lhuiller, F. Mila, M. 
Horvatic, and B. Ducot. 



P PUk 


P2a P4/,, J (K) e 


.81(1) .79(1) 


.95(2) .81(1) 94(9) .82(2) 


.89(1) .89(1) 


.94(2) .86(1) 85(6) .81(2) 


.96(1) .96(1) 


1.00(2) .94(1) 86(6) 0.80(1) 



TABLE I. Nominal occupation of Cr sites and occupations of Cr(12A;), Cr(2a), and 
Cr(4/t,j) sites from neutron refinement. J and e are the couplings extracted from the high-T 
mean-field analysis for the three samples. 
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FIGURES 
FIG. 1. Typical ^^Ga field sweep obtained at 40.454 MHz for p = .96. ^From right to left: 150 
K, 120 K, 100 K. The zero shift reference position is 3.116 Tesla. Horizontal arrow indicates shift 
direction upon cooling. Inset: Crystal structure of ideal SrCrgGaaOig. The thick dashed lines 
show the hyperfine coupling paths of Ga(4/) nucleus to frustrated Cr moments. Light grey circles 
are oxygens. Cr(12A;) are arranged to form a distorted kagome network and are coupled through 
Cr(2a). The full structure is Cr(4/j,j)-Cr(4/i,j)/kagome-Cr(2a)-kagome/Cr(4/^j)-Cr(4/i,j), etc. 

FIG. 2. K~^ versus T down to 15 K, for all samples. Minor second-order quadrupole corrections 
have been performed. The straight line extrapolation of K^^ = yields @nmr- 

FIG. 3. Calculated mean-field susceptibilities Xi2k and X2a of the Cr(12fc) and Cr(2a) moments 
versus T, with increasing Cr(2a)— Cr(12/i;) coupling. J is fixed to 85 K and pi2k = P2a = 1- 

FIG. 4. The (e, J, 5) set of values which reproduce the linear variation of K^^ with T of the 
p = .96 sample. For a given value of e (abscissa) there is an associated value S (upper panel) and 
J (bottom panel). 
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